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1. Introduction 

One of the primary goals of the CERN Large Hadron Collider (LHC) is a thorough inves- 
tigation of the mechanism of electroweak (EW) symmetry breaking and, more specifically, 
the discovery of one or more Higgs bosons and the determination of their properties jl], |||. 
In this context, vector-boson fusion (VBF) has emerged as particularly interesting class of 
processes. Higgs boson production in VBF, i.e. the EW reaction qq — » qqh°, where the 
Higgs decay products are detected in association with two tagging jets, offers a promising 
discovery channel || and, once its existence has been verified, will help to constrain the 
couplings of the Higgs bosons to gauge bosons and fermions j|. QCD corrections to the 
total cross section for single Higgs boson production via VBF have been computed using 
the structure function approach O. Distributions at NLO accuracy have recently become 
available through fully flexible parton-level Monte Carlo programs such as MCFM and 
VBFNLO @, 0, §. 

The two-Higgs doublet model (THDM) predicts the existence of two neutral CP-even 
Higgs bosons, one neutral CP-odd Higgs boson, and two charged Higgs bosons which 
have selfcouplings as well as couplings to gauge bosons and fermions || 10 1. Studies 



have shown that Higgs pair production at the LHC can serve as a probe of the Higgs 
potential [11, 12, 13]. Of interest recently has been the process pp — > h°h°jj — » bbbbjj 



via VBF in the context of the two higgs doublet model 14, 15, 16, [17], ffg. It was shown 



in Ref. [14] that in favorable THDM scenarios that it may be possible to extract the 
H° — > h°h° — > 46 resonance, thereby making the measurement of the trilinear h°h°H° 
coupling possible at the LHC. Assuming such a favorable scenario, the knowledge of QCD 
radiative corrections for this process will be needed in order to reduce the theoretical 
uncertainty on the total cross section and distributions. It is the aim of this paper to 
present the NLO QCD corrections for the process pp — > hPhPjj —* bbbbjj via VBF in the 
form of a fully flexible partonic Monte Carlo program within the type II-CP conserving 
two higgs doublet model. 
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In Section || we give an overview of the THDM and establish benchmark points used 
in our simulations. Section || lays out the details of the NLO calculation. Cross sections 
and distributions for the LHC are given in Section ||. Conclusions are given in Section ||. 



2. Two Higgs Doublet Model Parameters 



The THDM contains two SU(2) doublets, <3?i and <I>2i of weak hypercharge Y = 1. In the 
CP conserving THDM there is freedom in the choice of the Higgs boson-fermion couplings 
Q: type I, in which only one Higgs doublet couples to fermions; and type II, in which the 
neutral member of one Higgs doublet couples to up-type quarks and the neutral member of 
the other Higgs doublet couples to down- type quarks and leptons. Flavor changing neutral 
currents (FCNC) mediated by the Higgs bosons are automatically absent in both type I and 
type II THDMs |l9| . In this work, we only consider type II Higgs boson-fermion couplings. 
The most general THDM scalar potential which is invariant under SU(2)l <g> £7(l)y and 
conserves CP is given by |J , 
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with two real parameters, vi, V2 of mass dimension one and 6 real dimensionless parameters, 
Ai, . . . , A6- The minimum of the potential given by Eq. (|2.l| ) occurs at <I>j = (0, for 
(i = 1, 2). The physical spectrum of the Higgs sector of the CP-conserving THDM consists 
of two neutral CP even Higgs bosons, (h°, H°) , one neutral CP-odd Higgs boson, A , and 
two charged Higgs bosons, (H + , H~). The 8 parameters of the Higgs sector can also be 

taken as the vacuum expectation value, v = \J v \ + v \ = \J (\/2Gf)~ 1 , the masses of the 
Higgs bosons, m h o,m H o,m^o and m H ±, the mixing angles a and 0, and A5. 

The parameters of the Higgs potential, Eq. ( |2.1| ), can be restricted by imposing the- 
oretical requirements for the consistency of the model. We use, both, the requirement of 



vacuum stability [20] and perturbative unitarity |2l| for the tree-level coupling constants. 



The condition of vacuum stability is given by Eq. (2) of Ref. [22]. Perturbative unitarity 
requires that the magnitudes of all tree-level S-wave amplitudes for elastic scattering of 
longitudinally polarized gauge and Higgs bosons stay in the limit set by unitarity. Here, 



we consider the 14 neutral channels of Ref. [23]. p parameter constraints from electroweak 



precision data have also been considered |24j, [2q] . 

In this paper, we will consider two benchmark points tabulated in Table |l| which in fact 
satisfy the above mentioned requirements of vacuum stability and perturbative unitarity. 
Benchmark point Bl corresponds to a scenario in which the light Higgs boson does not 
couple to gauge bosons, i.e, cos(q — (3) = 1. In this scenario it will not be possible 
to produce a single light Higgs via VBF since the h°VV coupling is zero. Benchmark 
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point B2 is scenario in which the heavy Higgs H decouples from the gauge bosons, i.e., 
sin(a — (3) = 1. 





sin a 


tan (3 


A 5 




m h o 


m H o 


m H ± 


Bl 


0.832 


1.50 


-3.50 


295 GeV 


120 GeV 


300 GeV 


385 GeV 


B2 


-0.554 


1.50 





295 GeV 


120 GeV 


300 GeV 


385 GeV 



Table 1: THDM benchmark points. 



3. The NLO Calculation 

At leading order (LO), light Higgs pair production via VBF can, effectively, be viewed 
(see Fig. |l]a) as the elastic scattering of two (anti)quarks, mediated by i-channel W or Z 
exchange, with two light Higgs boson radiated off the weak-boson propagator. The "blobs" 
in Fig. [l] for the process qQ — * qQh°h° represent the vector boson scattering processes 
W+W~ -» h°h° (see Fig. §) and ZZ -» h°h° (see Fig. |), for charged current (CC) 
and neutral current (NC) processes, respectively. The generalization to crossed processes 
(q — > q and/or Q — > Q) is straightforward. In principal we should consider the double 
higgstrahlung process Vh h with V — > qq and the exchange of identical fermions in the 
initial or final state. However, in phase space regions with widely separated quarks jets of 
high invariant mass, the interference of these additional graphs is strongly suppressed by the 
large momentum transfer in the weak-boson propagators. Color suppression further makes 
these effects negligible. We, therefore, treat double higgstrahlung as a separate process 
and systematically neglect any identical particle effects as in the case of single Higgs boson 
production via VBF [^6|, g|. Further, gluon fusion h°h°jj production is treated as a separate 
process, since, the 0{a 2 a 2 s ) corrections are on the order of atobarns [27, p8f[ . 



The NLO calculation is performed in complete analogy to Ref. ||. The real emission 
graphs can be obtained by attaching the gluon to the quark lines of Fig. [l|a in all possible 
ways. Two distinct non-interfering color structures result: Feynman graphs with a single 
gluon attached to the upper quark line and Feynman graphs with single gluon attached 
to the lower quark line. Gluon initiated processes are obtained by crossing the final state 
gluon with the initial state (anti) quark. The result are graphs with i-channel and s -channel 
weak boson exchange. For consistency of the calculation we neglect the s-channel process 
gq — ► Vh h q, since, we have neglected double higgstrahlung at LO. 

All amplitudes are calculated numerically, using the helicity-amplitude formalism of 
Ref. fl2t|. Matrix elements for VBF proceses take the general form, 



M« Q ~ T^yJp?, (3.1) 
with Jfi and J® being the two quark currents shown in Fig. |l[ We have used the THDM 



implementation of FeynArts and FormCalc [30, |3l| to generate model predictions for the 
tensors T^ w and T^ v z for CC and NC vector boson fusion processes, respectively, for the 
case of off-shell vector bosons. We have introduced finite widths for h and H° in the 
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s-channel Higgs propagators. The subsequent decay of the light Higgs boson hP to bb is 
performed within the narrow-width approximation. 

Divergences arising from the real corrections are regulated in d = 4 — 2e spacetime 
dimensions by using the Catani- Seymour dipole subtraction method [32] in the dimensional 
reduction scheme [33]. Formulas for the subtraction terms and finite collinear pieces are 
identical to the ones for single Higgs production via VBF and are given in Ref. 0. 



V 
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h° 

~. h° 



<1 
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Figure 1: Feynman graphs contributing to qQ — * qQhPh at (a) tree-level and (b) including 
virtual corrections to the upper line. The "blobs" correspond to VVhPhP effective vertices which 



are represented by the tensors Ty V where V 



Z, W ± . 



4. Results for the LHC 

The goal of our calculation is a precise prediction of the LHC cross section for light Higgs 
boson pair production in VBF with two or more jets in the context of the THDM. In 
order to reconstruct jets from the final-state partons, the kx algorithm [53] as described in 



Ref. [35] is used, with resolution parameter D = 0.8. These jets are required to have 

p Ti >20GeV, \yj\<4,5. (4.1) 

Here yj denotes the rapidity of the (massive) jet momentum which is reconstructed as the 
four- vector sum of massless partons of pseudorapidity |r?| < 5. 

At LO, there are exactly two massless final state partons. The two hardest jets are 
identified as tagging jets, provided they pass the kx algorithm and the cuts described above. 
At NLO these jets may be composed of two partons (recombination effect) or three well- 
separated partons may be encountered, of which at least two satisfy the cuts of Eq. ((OJ) 
and would give rise to either two or three-jet events. As with LHC data, a choice needs to 
be made for selecting the tagging jets in such a multijet situation. Here the "py-method" 
is chosen. For a given event, the tagging jets are defined as the two jets with the highest 
transverse momentum with 

p£f > 30 GeV, \yf s \ < 4.5. (4.2) 

b -jets arising from decays of the two light Higgs bosons (h — > bb), are restricted by the 
following cuts: 

PTb > 30 GeV , \rj b \ < 2.5 , AR jb > 0.6 , AR bb > 0.7 (4.3) 
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Figure 2: Feynman graphs for the process W + W — ► h°h°. 



where ARjf, denotes the jet-6 separation in the rapidity-azimuthal angle plane. In addition, 
the 6-jets are required to fall between the two tagging jets in rapidity: 

l&n < y^max ■ (4-4) 

Backgrounds to VBF are significantly suppressed by requiring a large rapidity sepa- 
ration for the two tagging jets. Tagging jets are required to reside in opposite detector 
hemispheres with 

y tag 1 . y tag 2 < Q (4 5) 

and to have a large rapidity separation of 

Ay„ = |yf g 1 - yf* 2 | > 4 , (4.6) 
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Figure 3: Feynman graphs for the process Z°Z° — > 



sometimes called "rapidity gap cut". QCD backgrounds can be reduced by imposing a 
lower bound on the invariant mass of the tagging jets of 



(pf S 1 + Pf g 2 ) 2 > 600 GeV. 



(4.7) 



In all subsequent calculations we use the input parameters for defining Standard 
Model (SM) couplings as listed in Table ||. Other SM couplings are computed using LO 
electroweak relations. Cross sections are computed using CTEQ6M parton distributions 
||] with a s (M z ) = 0.118 for all NLO results and CTEQ6L1 parton distributions with 
a s (Mz) = 0.130 for all leading order cross sections. The running of the strong coupling is 
evaluated at two-loop order for all NLO results. In the following we use benchmark point 
Bl for NLO and LO Monte Carlo simulations. 

In Figures [|] and ||, we show the scale dependence of the total cross section within 
the cuts of Eqs. Q-Q for the 

process pp — > h°h°jj — > bbbbjj via VBF at the LHC. 
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Scale variations are shown for (a) LO results with (j,f = C/^o (black dotted line), (b) NLO 
results with /ip = £/^o and /j,r = /j,q (blue dot-dashed line), (c) NLO results with /jlr = £/xo 
and fip = Mo (green dashed line), and (d) LO results with fj,p = ^/iq (red solid line). In 
Figure |3], we choose a fixed reference scale, //o = 2m^. In Figure [| we choose the reference 
scale to be the virtuality of the exchanged t-channel vector boson, qy, shown in Figure [I]. 
For up = £// with 0.2 < £ < 5 the scale variation of the LO cross section is +19% to —14% 
for both choices of reference scale hq. While at NLO, the scale variations are —1.34% to 
-2.64% for fj, F = n R = 5 ±l {2m h ) and -5% to 0.5% for fj, F = fj, R = 5±V- Table | lists 
total cross sections at NLO for the benchmark scenarios of Table [lj For benchmark B2 
the cross section is below 1 femtobarn due to the fact that the heavy Higgs boson does not 
couple to gauge bosons. 

Our Monte Carlo program allows the analysis of arbitrary infrared and collinear safe 
distributions with NLO QCD accuracy. In order to assess the impact of the NLO corrections 
we compare LO and NLO results by plotting the dynamical K factor 



K(x) 



da^ LO (n R = ji F = no)/dx 
da%°((i R = fi F = fj, ) / dx 



(4.8) 



for our reference scale of hq. The stability of the results is represented via the scale 
dependence, given by the ratio of cross sections, 



R(x) 



dcr 3 ((i R = = CfJ-o)/dx 
da 3 (/ii 



= HF = Ho)/dx 

qy for NLO and LO distributions. 



(4.9) 



,p t T ag2 : 



Benchmark 


CTL0(lb) 


ONLo(fb) 


-KT-factor 


Bl 
B2 


86.87 
0.1864 


87.62 
0.1858 


1.0086 
0.9969 



We plot results for £ = 1/2 and 2 with // 

Figure ^ shows the maximum trans- 
verse momentum of the two tagging 
jets defined as Ppl^g = max(p^ s 1 
The left panel shows the distribution 
at LO (dashed green histogram) and 
at NLO (solid red histogram). The 
right panel shows phase space depen- 
dent K factor (solid green histogram) . 
Scale variations fiR = fip = £^0 for £ = 0.5 and 2 are shown for LO (dotted histograms) 
and NLO (dashed histograms). While at LO the scale variations are ±8%, they are reduced 
to ±1% at NLO. Figure || shows the rapidity separation of the two tagging jets. At NLO 
the tagging jet rapidity separation tends to increase which is typical behavior for VBF 
production processes ||. 



Table 3: Integrated LO cross sections, ulo, NLO cross 
sections, <jlo, in femtobarns (fb). The QCD iif-factors 
are defined as K = <jnlo I vlo- The renormalization 
and factorization scale has been set to /in = fip — qy- 



M z 


M w 


Gp 


91.188 GeV 


80.416 GeV 


1.16639 x 10" 5 /GeV 2 



Table 2: SM input parameters 
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Figure shows the distribution in the maximum transverse momentum of the four b- 
jets resulting from the decay of the light Higgs bosons defined as p™^ x = max(px,6l> • • • ,PT,bi) 
The distribution is peaked in the vicinity of 100 GeV because we have chosen the light Higgs 
boson mass to be rrih = 120 GeV. A lower light Higgs mass would result in a softer 6— jet. 
Figure |§] shows the distribution in the four 6-jet invariant mass, = \J {jpn + • • • + Pbi) 2 - 
A peak occurs a around the mass of the heavy Higgs boson, m# = 300 GeV. In both dis- 
tributions the scale dependence is reduced with relatively flat phase space dependent K 
factors. 
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Figure 4: Scale dependence of the total cross section at LO and NLO within the cuts of Eqs. (4.1)- 
(|4.7|) for VBF h°h°jj production at the LHC. The factorization scale fip and the renormalization 
scale [ir are taken as multiples, £/io, of the fixed reference scale fio — 2m,h- The NLO curves are 
for — /if = Cmo (solid red line), /if = Mo and [ir = £//o (dashed green line), and \xp = £/^o and 
Mfl — Cmo (dot-dashed blue line ). The dotted black curve shows the scale dependence of the LO 
cross section for [ip = £/io- 



In Figure |To|, we show the normalized azimuthal angle correlation of the tagging jets, 
A(f)jj, for benchmark points Bl and B2. Here A^-j = (j)j + — 4>j_ with j+(j-) being 
the "toward" ( "away" ) jet as defined in Ref. |37]]. For benchmark Bl, the normalized 
distribution in Ac/>jj has the characteristic shape for single H° production followed by the 
decay, H° — > h°h° where the tensor structure of the HVV vertex is g^ u |37], ^8|. However, 
for benchmark B2 in which the H° is decoupled from the gauge bosons the shape of the 
A<f)jj distribution develops a peak at A(j)jj = degrees as opposed to a dip. 1 This is due 
to the interference of s-channel and t-channel VV — ► h°h° graphs. 



1 Such features have been pointed out in Ref. |39| for slepton pair production via VBF. 
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Figure 5: The caption is the same as in Fig. except /iq = qv- 
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Figure 6: Maximum transverse momentum of the two tagging jets. In the left panel, der/dp™^ is 
shown at LO (dashed green) and NLO (solid red) for ^ = Qv ■ The right-hand panel depicts the K 
factor (solid green) and scale variations of LO (dotted) and NLO (dashed) results for fj,R = \if = £/ito 
with f = 1/2 and 2. 



5. Conclusions 

We have computed the next-to-leading order QCD corrections for light Higgs pair pro- 
duction via vector boson fusion at the LHC within the type II CP conserving two-higgs 
doublet model. We have included the subsequent decay of the light Higgs to the bb final 
state. Our NLO calculation takes the form of a fully flexible partonic Monte Carlo program 
allowing arbitrary phase space cuts. We have shown that scale variations for total cross 
sections and distributions are reduced at NLO. QCD if -factors are modest. These results 
are consistent with those of Standard Model Higgs production via vector boson fusion Q. 
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Figure 7: Maximum transverse momentum of the four b jets. In the left panel, da/dp™ 3 ^ is shown 
at LO (dashed green) and NLO (solid red) for /xq = qv- The right-hand panel depicts the K factor 
(solid green) and scale variations of LO (dotted) and NLO (dashed) results for hr = (if = £/ito 
with f = 1/2 and 2. 
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Figure 8: The four 6— jet invariant mass, m\ b = (pb,i + Pb,2 + Pb,3 + Pm) 2 ■ m the left panel, 
da/dm4b is shown at LO (dashed green) and NLO (solid red) for fi = Qv- The right-hand panel 
depicts the K factor (solid green) and scale variations of LO (dotted) and NLO (dashed) results 
for = \ip = £/io with £ = 1/2 and 2. 

We also note the sensitivity of the azimuthal angle correlation of the two tagging jets, 
da/dA(j)jj, to the tensorial structure of VV — ► h h scattering amplitudes. For THDM 
scenarios in which the heavy Higgs coupling to electroweak gauge bosons is highly sup- 
pressed, the Acpjj distribution is peaked at A(f>jj = degrees while for scenarios in heavy 
Higgs couples to electroweak gauge bosons there is a dip at A0jj = degrees. 
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Figure 9: Rapidity separation within the cuts of Eqs. (4.1)-(4.7). In the left panel, da/dyjj is 
shown at LO (dashed green) and NLO (solid red) for fXo = qv- The right-hand panel depicts the K 
factor (solid green) and scale variations of LO (dotted) and NLO (dashed) results for ^lr = ^f= £a*0 
with f = 1/2 and 2. 
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Figure 10: The normalized distibution in azimuthal angle correlation of the two tagging jets, 
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, for benchmark scenarios Bl (red) and B2 (blue) at NLO (solid) and LO (dashed). 
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